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Electron sheaths form near the surface of objects biased more positive than the plasma potential, 
such as a Langmuir probe collecting electron saturation current. Generally, the formation of electron 
sheaths requires that the electron-collecting area be sufficiently smaller (^2.3me/M times) than 
the ion-collecting area. They are commonly thought to be local phenomena that collect the random 
thermal electron current, but do not otherwise perturb a plasma. Here, using experiments on an 
electrode embedded in a wall in a helium discharge, particle-in-cell simulations, and theory it is 
shown that under low temperature plasma conditions (Te ^ D) electron sheaths are far from local. 
Instead, a long presheath region (27 mm, approximately an electron’s mean free path) extends into 
the plasma where electrons are accelerated via a pressure gradient to a flow speed exceeding the 
electron thermal speed at the sheath edge. This fast flow is found to excite instabilities, causing 
strong fluctuations near the sheath edge. 


INTRODUCTION 

A sheath is the space charge region found near the 
physical boundaries of most plasmas. The vast major¬ 
ity of sheaths are ion rich because this is what natu¬ 
rally forms as highly mobile electrons charge a surface 
negatively. Comparatively little is known about electrons 
sheaths, although they are routinely produced when ob¬ 
jects are biased positive with respect to the plasma po¬ 
tential [T]. 

The most common situation is the electron satura¬ 
tion region region of a Langmuir probe sweep, but they 
arise in many other situations including negative ion 
sources [2] and electron sources [3], positive electrodes 
employed for blob control [4] , particle circulation in dusty 
plasmas [5], and turbulence-induced particle fluxes [6]. 
Electron sheaths are also common in several other sit¬ 
uations, including: near highly emitting surfaces 0 , in 
microdischarges [8], during the high potential phase of 
the rf cycle in processing discharges [9], around electro¬ 
dynamic tethers HO], the lunar photosheath m, around 
wire arrays used for electron temperature control [12], 
and in scrape off layer control m- 

In this Letter, results of experiments, particle-in-cell 
(PIC) simulations and theory are provided showing that 
electron sheaths form a long electron presheath extending 
well into the plasma. Eurthermore, electrons are acceler¬ 
ated to high velocities in this region, obtaining a distribu¬ 
tion that is flow-shifted to a speed exceeding the electron 
thermal velocity at the sheath edge. This may be consid¬ 
ered an electron sheath analog of the Bohm criterion [T4| . 
The fast differential streaming between electrons and ions 
is found to excite streaming instabilities that give rise to 
strong fluctuations of the boundary layer region. 

Electron presheaths are found to differ in their essen¬ 
tial properties from ion presheaths. In particular: (1) the 
differential potential (Acf) is much smaller, nominally by 
a factor of Ti/T^, (2) it is much longer in extent, nomi¬ 


nally by a factor of ^/rui/rrie, (3) electrons are acceler¬ 
ated by a pressure gradient, in contrast to direct electric 
field acceleration of ions in an ion presheath, and (4) 
the differential streaming excites instabilities and strong 
fluctuations. These results promise new insights into the 
applications mentioned above. Advances in the funda¬ 
mental physics of electron sheaths may also lead to new 
applications. 


EXPERIMENTS 

Experiments were conducted in a cylindrical vacuum 
chamber in which an electrode with a diameter of 19 
mm was embedded into one face. The walls and faces 
of the vacuum chamber were maintained at ground po¬ 
tential. The area ratio of the chamber (1.86 x 10^ mm^) 
to the electrode (1.14 x 10^ mm^) was chosen to ensure 
the formation of an electron sheath above the auxiliary 
electrode when biased above the plasma potential m- 
A plasma was generated in 20 mTorr of helium with a 
barium-impregnated tungsten thermionic emitter (Heat¬ 
wave Labs model 101117) located 10 cm from the elec¬ 
trode. The emitter was operated at a temperature of 
approximately 1.1 x 10^°C. A cross section of the exper¬ 
imental setup can be seen in Eig. 

The discharge current was held constant at 300 mA 
which resulted in a discharge voltage which varied from 
49 to 54 V over the course of the measurements. The 
auxiliary electrode was biased to -50 V, 0 V and 15 V, 
forming an ion sheath, a weak ion sheath, and electron 
sheath respectively. An emissive probe was operated at 
a height of 35 mm and at a radial position of 70 mm and 
was used to measure the plasma potential via the float¬ 
ing point technique. The plasma potential was 6.2 V, 5.3 
V, and 6.6 V for the three aforementioned cases. Two- 
dimensional maps of the electron densities above the elec¬ 
trode were generated using the laser-collisional induced 
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FIG. 1. Sketch of the experimental setup used to measure the 
electron densities above an embedded electrode. 
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FIG. 2. Measured electron density prohles above the electrode 
along the center axis at -50 V (solid line), 0 V (dotted line), 
and 15 V (dashed line). The profiles are normalized by their 
values at 35 mm: 3.5 x 10^ cm“^ at -50 V, 3.3 x 10^ cm“^ at 
0 V, and 3.7 x 10^ cm“^ at 15 V. Also included is a gray line 
representing the constant density gradient found at z > 30 
mm. 


fluorescence (lcif) diagnostic [16]. High energy electrons 
could potentially cause secondary electron emission, how¬ 
ever the uncollided electron density (high energy tail) is 
estimated to be two orders of magnitude below the mea¬ 
sured densities. Given a secondary emission coefficient of 
7 ~ 0.1, the density of electrons from secondary emission 
should be about three orders of magnitudes below the 
measured densities. 

Figure shows the axial density profiles from LCiF 
measurements (relative uncertainty estimated to be 
7% Hlj), normalized by their density at 35 mm, for the 
ion sheath (-50 V, solid line), a weak ion sheath (0 V, 
dotted line), and the electron sheath (15 V, dashed line). 
In all three cases the same linear density gradient is ob¬ 
served far from the electrode (> 30 mm). 

This gradient is likely the result of the nonuniform gen¬ 
eration of plasma in the volume. Attenuation of the pri¬ 
mary electrons and their falloff with distance from the 
small source will lead to a locally peaked density distri¬ 
bution. This non-uniformity will lead to diffusion to the 
walls of the chamber and an associated ambipolar field. 
Based on Langmuir probe measurements of the electron 
temperature and observed density gradients this ambipo¬ 
lar field will be of the order 0.6 V/cm and will tend to 
inhibit the transport of electrons toward the embedded 


electrode. 

At -50 V and 0 V, the electron depletion region asso¬ 
ciated with an ion sheath can clearly be seen extending 
from 0-4 mm and 0-3 mm respectively. The presheath 
is estimated to be the region at which the electron den¬ 
sity gradient deviates from the roughly constant slope 
(> 30 mm). The estimated location of the presheath 
edge occurs near 12 mm. This can be compared to the 
calculated ion-neutral mean free path which is expected 
to determine the presheath length scale [T8|. Estimating 
ion energies of 0.1 eV and using the isotropic scattering 
cross sections of Phelps [19] (a = 2.4 x 10“^^ m^), we 
find a mean free path of 6.5 mm which is on the order of 
the estimated location of the presheath edge. 

The electron sheath (15 V) also features a region of 
electron depletion from about 0-2 mm from the electrode. 
This is caused by the acceleration of the electron fluid 
and its consequent rarefaction. The salient feature of the 
electron sheath case is that its density gradient is essen¬ 
tially the same as the ion sheath’s far from the electrode 
(> 27 mm), but closer to the electrode it bears a dis¬ 
tinctly different gradient. Assuming that this deviation 
is caused by collisional effects similar to the ion sheath 
case, a calculation can be made for electrons based on 
their momentum transfer cross section. Assuming an en¬ 
ergy of 4.0 eV (a = 6.6 x 10“^^ m^ per Phelps [E]), the 
electron mean free path is found to be 24 mm. 

This suggests that the positively biased electrode is af¬ 
fecting the plasma far from the electron sheath region, on 
the order of an electron mean free path, and suggests the 
presence of an electron presheath. The Debye length at 
the sheath edge is estimated to be 0.38 mm based on the 
LCIF measurements and Langmuir probe measurements 
of the electron temperature. Thus this presheath extends 
approximately 70 Debye lengths from the boundary. In 
both the case of the electron presheath and ion presheath, 
the length scales are significantly smaller than the length 
scale of the system. 


SIMULATIONS 

A complementary analysis was carried out using Aleph, 
a PIC-DSMC code developed at Sandia National Lab¬ 
oratories. Aleph is an electrostatic code intended for 
massively paralllel (over 10k cores) plasma simulations. 
Fields are solved using finite-element method libraries 
from the Trilinos project [20]. Particles are advanced 
using the velocity Verlet algorithm as described by Spre- 
iter and Walter m- Other recent uses of Aleph include 
the study of vacuum arc discharges [22] and the onset of 
plasma potential locking [23] . 

The simulation was conducted in a rectangular Cartes- 
tian domain, 75 mm by 50 mm, with an unstructured 
mesh formed of triangles. This is the same geometry as 
that used in m and depicted in Figure]^ The left-hand 
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FIG. 3. Representation of the simulation geometry with over¬ 
laid electric field potential contours from electron sheath sim¬ 
ulations. The domain is 75 mm by 50 mm with an axis of 
symmetry on the left hand side. 


side of the domain specularly reflects all particles and 
possesses a Neumann boundary condition of dVjdx = 0, 
thus representing an axis of symmetry. Aside from the 
electrode shown in Figure the outer boundaries are 
held at R = 0 and out flux all incoming charged particles. 
The characteristic size of the triangles used to mesh the 
domain was 233 /im which resolved the Debye length of 
approximately 540 jam. All simulations used a timestep 
of 25 ps, which is sufficient to meet CFL requirements 
for all particles observed in the simulation and resolves 
the plasma frequency of 0.28 GHz. 

The simulation domain maintains an area ratio of wall 
to electrode similar to that in the experiment, chosen to 
assure the formation of an electron sheath [HI [25] . As 
the formation of the electron sheath only depends on the 
area ratio and not the absolute size of the electrode, ex¬ 
periment and simulation should be largely comparable. 
The simulated electrode was biased to -50 V and 15 V to 
match the measured ion and electron sheath. Quasineu¬ 
tral plasma (Tg = 4.0 eV, = 1000 K) was added to 
the simulation domain at a constant rate in a rectangu¬ 
lar area located 37.5 mm from the anode. Simulations 
were run for 30 /is at which point they were found to 
be in equilibrium based on field energy and total parti¬ 
cle number. Field and particle properties were averaged 
over an additional 20 /is in order to minimize statistical 
fluctuation in quantities of interest. 

The experimental density measurements are compared 
to the equivalent simulations in the plots on the left side 
of Fig. Overlaid on the simulated electron densities 
are arrows showing the electron flux vectors scaled to 
the same value for both cases. The horizontal axes have 
been normalized by the electrode radii, in order to pro¬ 
vide a more suitable comparison. The right side of Fig.|^ 
presents maps of the charge density normalized by the 
density of the collected species. Inset in the charge den¬ 
sity maps are power spectra of the sheaths’ positions over 
time. 

In the case of the ion sheath, a large region of elec¬ 
tron depletion is visible above the face of the electrode. 


The size and shape of this region is largely consistent 
between simulation and experiment, with the small dif¬ 
ferences likely ascribable to a discrepancy in the electron 
densities at the sheath edge. An ion sheath is also ob¬ 
served above the grounded wall. Electron current is small 
throughout the simulated domain and is consistent with 
the fact that ion sheaths tend to confine electrons. The 
charge density shows the formation of a stable sheath. 

The electron sheath simulation also features a region of 
electron depletion near the face of the electrode, resulting 
from acceleration of the electron fluid. Though the width 
of the sheath in the density maps from simulation appears 
qualitatively larger in the density maps, the sheath edge 
seen in the charge densities oscillates about 2 mm consis¬ 
tent with the experimental electron density profile seen 
in Fig. Further differences are likely attributable to a 
combination of potential causes: lower electron densities 
in the simulation and the absence of electron-neutral col¬ 
lisions. The former would result in larger Debye lengths 
and subsequently larger sheaths. The latter would also 
tend to increase the sheath width as seen in equation (6) 
of [26] . 

The electron sheath simulation also features a substan¬ 
tial degree of electron current directed toward the elec¬ 
trode from at least as far away as 30 mm. The plasma 
properties in this region (rig = 4.5 x 10^ cm“^, Tg = 2.4 
eV) give a Debye length of 0.54 mm, 60 times smaller 
than the extent of the directed flow. While it is well 
known that ion presheaths can extend throughout the 
entire plasma for large mean free paths m, we will 
show that (for the same collision process) the electron 
presheath is nominally longer. The extent of the flow 
is consistent with where the density profile from Fig. 
begins to change slope, indicating the presheath. The 
electron presheath is clearly a significant perturbation to 
the plasma. Furthermore, as opposed to the ion sheath, 
the electron sheath edge exhibits significant fluctuations 
in its position. 

The proximity of the grounded wall to the auxiliary 
electrode leads to a funnel-like structure in the electron 
flow with a notable convergence. We note that there is 
the possibility for some overlap of the electron presheath 
with the ion presheath from the abutting walls. In this 
case, the radial electric field of the ion presheath will can¬ 
cel as a result of the symmetry of the system. The axial 
electric field will tend to accelerate electrons away from, 
rather than toward the anode, therefore the flows ob¬ 
served in Fig.j^are not the result of adjacent ion sheaths. 
Finally, we note that the ion presheath length scale is es¬ 
timated to be 6.5 mm, significantly less than the observed 
electron presheath size and the device dimensions. 
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FIG. 4. Experimental and simulated electron density maps (left) and simulated on-axis charge densities over time (right). 
The ion sheath results comprise the upper plots, and the electron sheath results comprise the lower ones. Overlaid atop the 
simulated electron densities are the electron particle flux, the magnitudes of which range from 0.1-1.2x10^° m“^s“^ for the 
electron sheath, and below 5 x 10^ m“^s“^ for the ion sheath. The power spectral densities of the sheath edge position are 
inset in the charge density plots. 


THEORY 

Traditional Langmuir probe analysis assumes that a 
probe in electron saturation collects the random ther¬ 
mal flux of electrons incident on the electron sheath. An 
implication of this local picture is that the electron ve¬ 
locity distribution function (evdf) at the edge of the 
electron sheath would be a half-Maxwellian with no flow 
shift (but with a flow moment) [27H29] . The random flux 
of electrons flows into the sheath from the quasineutral 
plasma, and all electrons reaching the sheath are lost to 
the boundary. However, our results show that the pres¬ 
ence of an electron presheath leads to a vastly different 
picture. 

In particular, simulations show that the presheath is 
found to introduce a substantial flow-shift in the electron 
distribution, approaching the electron thermal speed by 
the sheath edge; see Fig. Contrary to the conven¬ 
tional picture of a highly-kinetic truncated distribution 
function, this figure shows that it is in fact well repre¬ 
sented by a flowing Maxwellian. That the distribution 
is well-described by a Maxwellian despite the absence 
of electron-neutral collisions is notable. Confirmation of 
this flow shift in experiment remains an open challenge. 
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FIG. 5. The electron velocity distribution functions normal 
to the electrode, along x = 0, at several locations above the 
electrode in the simulations: the injection region (a), the 
sheath edge (o), and inside the sheath (o). 


A significant number of electrons flow out of the electron 
sheath (negative velocities are electrode-directed) despite 
the absence of an explicit collision algorithm in the sim¬ 
ulations. The sheath is defined as where quasineutrality 
is sufficiently violated, or [n^—n^jn^ = e where we have 
chosen e = 0.3 in order to avoid stochastic density vari¬ 
ations in the plasma. This approach places the sheath 
edge at z = 2.0 mm. More is said of this choice for e 
below, in view of the strong density fluctuations in this 
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region. 

While PIC simulations can result in numerical ther- 
malization of nonequilibrium distributions, this is not ex¬ 
pected to be a factor in the present results. Estimates 
based on the work of Montgomery and Nielsen [30] sug¬ 
gest that the thermalization time for the present system 
is of the order 2.6 /is. This exceeds the 1.1 /is required 
for the electron fluid to transit from the source region 
to the electrode. Additionally, recent simulation results 
of the ion to electron sheath transition [31] using com¬ 
parable simulation parameters possess evdfs exhibiting 
both a flow shift and a loss cone due to the presence of 
the boundary. Both factors are found to affect the flow 
speed, with the flow shift being at least as important as 
the loss cone. 

A fluid analysis is used to interpret aspects of the ex¬ 
periments and simulations and to better understand the 
origin of the observed flow shift. Consider the fluid mo¬ 
mentum balance 


Ue 


dUe 

dz 


eE 

me 


k^Te dfle 
mefle dz 


- Ue{yc + T^s) 


( 1 ) 


where e is the elementary charge, E is the electric field, 
Vc is the collision frequency, and Uq represents the source 
(ionization) frequency. The terms on the right-hand side 
represent the forces due to the electric field, pressure gra¬ 
dient, and collisions respectively. The largest of these 
terms is the pressure gradient term. 

The ions are estimated to vary according to the Boltz¬ 
mann density relation. 


eE 1 dn 
ksTi n dz 


( 2 ) 


This assumption depends on a number of factors which 
may not necessarily be met for all electron sheaths. 
Multi-dimensional effects, flow, and collisions may all re¬ 
duce the applicability of this relation. Indeed, in the sim¬ 
ulations conducted, this relation is not strictly applicable 
as the ions experience substantial collisional and inertial 
forces as seen in Figure The figure depicts the iner¬ 
tial (solid), pressure (dotted), electric field (dash-dotted), 
and collisional (dotted) fluid terms for the ions along the 
axis of symmetry, neglecting perpendicular components. 
These cases require significantly more complex analysis; 
an example of the treatment of the multi-dimensional 
case can be found in [26]. Therefore, this relation is as¬ 
sumed here strictly as a means of developing an initial 
estimate of electron sheath and presheath properties in 
the ideal case of minimal ion-neutral collisions and neg¬ 
ligible flow. We note that no such assumption is made in 
the simulations. 

Substituting this into Eq. 0 shows that the pressure 
gradient term is Te/Ti larger than the electric field. This 
is consistent with previous emissive probe measurements 
in a discharge of lower density (10^ cm“^ compared to 10^ 
cm“^) which showed essentially no field past 4 mm from 



FIG. 6. Calculation of the ion fluid terms along the axis of 
symmetry. Two-dimensional effects are excluded. The iner¬ 
tial term (RHS Eq. [^ is the solid line, the electric field term 
(first term on LHS) is the dash-dotted line, the pressure term 
(second term on LHS) is the dotted line, and collisions (third 
term on LHS) are shown by the dashed line. 


the electrode [25]. Although the potential gradient in 
the presheath is small (characterized by T^), the resulting 
pressure gradient drives a strong electron flow due to the 
density gradient that results. This contrasts with the 
situation found in ion presheaths, where the electric field 
term exceeds the pressure gradient term by Te/T^, and 
the presheath potential drop is of the order of Tg (rather 
than Ti). In this case, ions are accelerated ballistically by 
the electric field to a speed exceeding the sound speed at 
the sheath edge. The importance of the pressure gradient 
term with respect to the electric field is confirmed in 
the simulations. Figure [^a) plots the pressure gradient 
(solid line) and electric field (dashed line) acceleration 
terms from Eq. 0 on a logarithmic scale, as a function 
of distance from the electrode. 

The simulations indicate that the pressure gradient is 
the dominant acceleration mechanism from 2-5 mm and 
z > 11 mm. The drop in pressure gradient between 5 
mm and 11 mm coincides with a plateau in the simu¬ 
lated electron density near the sheath edge (similar to 
that seen in Fig. and a plateau in electron tempera¬ 
ture (not shown). Several factors may contribute to this 
plateau including a stagnation of ions as they approach 
the sheath potential barrier or a change in the conver¬ 
gence of the electron fluid. Past 11 mm, the pressure 
gradient continues to dominate the acceleration of the 
electron fluid up to the sheath edge. 

The degree to which the electron fluid is accelerated 
can be calculated via the electron continuity equation 


d 

dz 


(neUe) = VgUe, 


(3) 


and a common sheath criterion, which identifies the 
sheath edge as the location where ^qdn/dz < 0 
[18], where q is the charge of the species. Dropping 
the collision terms, Eqs. 0 and 0 yield dnejdz = 
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FIG. 7. (a) Pressure gradient (solid line) and the electric 

field (dotted line) accelerations obtained from simulations. 
The solid gray lines indicate the locations where the evdfs 
in Fig. were obtained, with the middle line indicating the 
location of the sheath edge, (b) Normalized charge density 
(dashed line), the electron fluid velocity (solid line), and po¬ 
tential (dash-dotted line) in the direction normal to the elec¬ 
trode surface, with respect to the distance from the electrode 
along X = 0.0. 


eneEf{meul — ksTe). From Eq. ([^, and assuming 
quasineutrality, the sheath criterion provides an electron- 
sheath analog of the Bohm criterion, 




lkB{Te^Ti) _ 

- = Vb- 


m. 


(4) 


A similar derivation of the electron fluid speed at the 
sheath edge was previously obtained by Loizu, Dominski, 
Ricci, and Theiler [32] This criterion demands a region 
of electron acceleration outside of the sheath region. 

Figure nib), shows the charge density (dashed line), 
the electron fluid velocity (solid line), and the local po¬ 
tential (dash-dotted line). The electron fluid is found 
to reach a velocity of 0.92'r’T by the sheath edge (as 
previously defined), in fair agreement with Eq. Q. A 
number of factors may contribute to the remaining dis¬ 
crepancy including the ambiguity in the definition of a 
precise sheath edge location [24| , the use of a planar one¬ 
dimensional theory in describing a converging flow, the 
presence of adjacent sheaths, and the substantial fluctu¬ 
ations observed in the simulations near the sheath edge. 

The fast electron flow creates a large differential 
streaming between electrons and ions that is expected 
to lead to electron-ion streaming instabilities [33] in the 
electron presheath and sheath. Indeed, Eig.j^shows sub¬ 
stantial fluctuations of the electron sheath edge, but not 
in the case of the ion sheath simulation. The frequency 
of these fluctuations is observed to peak around 0.8 MHz 
comparable to the most unstable mode of ion-acoustic in¬ 
stability (ion plasma frequency) which would be 1.4 MHz 


for = 1.7 X 10^ cm“^ (z = 2.0 mm). A more in-depth 
investigation of the ion density fluctuations shows good 
agrement with the ion-acoustic dispersion relation [26] . 
This observation is consistent with previous observations 
of ion acoustic waves excited by positive probes [34] and 
current fluctuations in measurements using a segmented 
electrode [25]. It should also be noted that the low fre¬ 
quency of these fluctuations suggest that they are dis¬ 
tinct from previous instabilities observed near electron¬ 
collecting interfaces [35H37| . 

In addition, it is observed that the ion flow obtains a 
radial component near the sheath edge, which may con¬ 
tribute to ion-ion two-stream instabilities. The presence 
of these strong fluctuations blurs the sheath edge loca¬ 
tion, as shown in Eig. In a steady-state sense, a tran¬ 
sition region is formed between presheath and sheath as 
a result of the fluctuations. The previous estimate of 
the sheath edge (at z = 2.0 mm) was based on where the 
time-average charge density is reduced by 30% (if instead 
an estimate based on the Child-Langmuir sheath thick¬ 
ness is used, Eq. (13) of [29], an estimate of z = 1.5 mm 
is obtained [38]). 

The length scale of an ion presheath is typically de¬ 
termined by collisional processes, and is estimated as 
the ratio of the ion flow speed to the collision fre¬ 
quency Cs/y for a specific process such as ionization. 
The electron presheath length would be estimated to be 
/e vtI^^ This implies that, for the same collision pro¬ 
cess, the electron presheath is much longer than the ion 
presheath, by a factor of 

_ /^e T ^2 ^ I /^\ 

li Cg V ^7/0 Tq y m/0 

As an example, le/h ~ 270 in argon or le/h ~ 85 in he¬ 
lium. However, this ratio is only applicable in the case 
of the same collision process. In the present case, it is 
believed that the ion and electron presheath lengths are 
likely governed by different processes, namely ion-neutral 
and electron-neutral collisions. This suggests le/U ~ 3 
based on calculations of the mean free paths. We also 
note that in typical low temperature plasma experiments, 
these lengths may be constrained by the dimensions of 
the plasma rather than collisions. Eigures [^ and sup¬ 
port the suggestion that electron presheaths are much 
longer than ion presheaths. This is a very different pic¬ 
ture than the assumption of a local phenomenon that is 
found in Langmuir probe theory. 


CONCLUSIONS 

These novel properties of the electron sheath are sur¬ 
prising both because of how they differ from ion sheaths 
and because of their influence on the bulk plasma. The 
perturbations in electron density and flow caused by what 
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would otherwise be considered a small electrode suggests 
that conventional models of electron sheaths need to be 
revisited. These fundamental physics results may also 
lead to useful new applications. 
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